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Temporal quark correlation functions are analyzed in quenched lattice QCD for two values of 
temperature above the critical temperature (T c ) for deconfinement, T = 1.5T C and 3T C . A two-pole 
ansatz for the quark spectral function is used to determine the bare quark mass and the momentum 
dependence of excitation spectra on large lattices of size up to 128 3 x 16. The dependence of 
the quark correlator on these parameters as well as the finite volume dependence of the excitation 
energies are analyzed in detail in order to examine the reliability of our analysis. Our results suggest 
the existence of quasi-particle peaks in the quark spectrum. We furthermore find evidence that the 
dispersion relation of the plasmino mode has a minimum at non-zero momentum even in the non- 
perturbative region near T c . We also elaborate on the enhancement of the quark correlator near the 
chiral limit which is observed at T = 1.5T C on about half of the gauge configurations. We attribute 
this to the presence of near zero-modes of the fermion matrix that are associated with non-trivial 
topology of the gauge configurations. 
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I. INTRODUCTION 

At asymptotically high temperature (T) strongly inter- 
acting matter, described by Quantum Chromodynamics 
(QCD), behaves like an almost free gas of quarks and 
gluons. Excitations of the medium in this temperature 
range are well described by dispersion relations of free 
quarks and gluons. As T is lowered, the temperature- 
dependent running coupling constant, g, increases and 
the excitation spectra of the elementary quark and gluon 
fields are modified. Eventually they deviate substantially 
from those of free particles. In particular at those low 
temperatures close but above the QCD phase transition, 
T c , at which experiments at the Relativistic Heavy Ion 
Collider (RHIC) Q found evidence for the existence of 
strong correlations in hot and dense matter, it is of inter- 
est to understand the fate of the elementary excitations, 
their quasi-particle properties and the structure of dis- 
persion relations. 

^From leading-order rcsummed perturbation theory it 
is known that collective excitations of the quark and 
gluon fields acquire mass gaps (thermal masses) and de- 
cay rates proportional to gT and g 2 T, respectively 0- 
Since at leading order the decay rates in units of T para- 
metrically grow faster than the thermal masses as T is 
lowered, it is naively expected that quasi-particle modes 
of these fields cease to exist at low temperatures at which 
g is large, even in the deconfined phase. On the other 
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hand, the experimentally observed quark number scaling 
of elliptic flow suggests the existence of quasi-particlcs 
having quark quantum numbers at the early stage of fire- 
balls created in heavy ion collisions Q . Also the behavior 
of quark number susceptibilities and higher-order cumu- 
lants of quark number fluctuations, calculated in lattice 
QCD, suggests the existence of quasi-particles even near 
T c 01. 

One of the striking features of the quasi-particle excita- 
tions of quarks, found in leading order perturbation the- 
ory, is that the in-medium quark dispersion relation splits 
into two branches, the normal and plasmino modes 
Moreover, the latter has a non-trivial minimum at non- 
vanishing momentum of order gT [f| . Such a spectrum of 
fermions is, in fact, quite common. In the weak coupling 
and high T limit it is realized not only for QCD, but also 
in a wide class of models where the fermion couples to 
a boson. Away from this limiting situation the fermion 
spectrum has been investigated in various theoretical set- 
tings 043- An interesting observation made in these 
studies is that the fermion spectrum exhibits a variety of 
structures depending on parameters of the system. For 
example, it has been pointed out that the fermion spec- 
trum has a multi-peak structure when the mass of bosons 
that couple to the fermion is comparable with T Q: In 
addition to the normal and plasmino peaks having ther- 
mal masses, there appears an additional, third dispersion 
which crosses the origin of the energy-momentum plane. 
The existence of such a multi-peak structure is also sug- 
gested in Refs. us, mm. 

It is thus interesting to inves- 
tigate the fate of this third branch in the quark spectrum 
in the non-perturbative region near but above T c . 

In order to gain insight into the quasi-particle nature 
of quarks in the non-perturbative region, analyses of the 
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spectral properties of quar ks have recently been carried 
out on the lattice QCD jl7H19j . Lattice calculations, per- 
formed in the quenched approximation of QCD in Lan- 
dau gauge [HEl for 1.25 < T/T c < 3.0, indicated that 
a plasmino mode with a distinctively different dispersion 
relation contributes to the quark spectral function in ad- 
dition to the normal mode even at temperatures close to 
T c . In the chiral limit these modes have identical quasi- 
particlc masses that are, in the range of temperature an- 
alyzed, approximately proportional to T. It is notable 
that similar results have been reported in the analysis 
of the quark spectrum with Schwinger-Dyson equations 

The lattice QCD calculations presented in Refs. [l7], 
HH, however, showed that the quasi-particle masses and 
dispersion relations of quarks are quite sensitive to fi- 
nite volume (V) effects. The quark thermal mass drops 
as the aspect ratio V X ^ 3 T = N a /N T increases, and it 
does not converge at the value on the largest lattice in 
Refs. [13, witri n <t/N t = 4, while the lattice spacing 
a used in these studies was found to be fine enough to 
suppress lattice discretization effects. The origin of such 
a strong V dependence is well understood in terms of dis- 
cretization effects in momentum space. For N a /N T = 4 
the lowest non-zero momentum is p m i n = 2tt(N t /N (T )T ~ 
1.6T for periodic boundary condition along the spatial 
direction, which is still significantly larger than T. On 
lattices with such an aspect ratio effects of thermally- 
excited particles with momenta of order T, which are re- 
sponsible for the emergence of the thermal mass in terms 
of perturbation theory 0, Ho| , will thus not be properly 
incorporated. 

In the present study, we analyze the quark correla- 
tion function on lattices with much larger aspect ratio, 
N a /N T = 8. With this aspect ratio, the lowest non- 
zero momentum, p mm ~ 0.79T, is lower than T, and 
effects of low momentum modes are better incorporated 
in numerical calculations. To extract the spectral func- 
tion from Euclidean lattice correlators, we use the two- 
pole ansatz as in the previous studies [Tt], E3 ■ We show 
that this ansatz reproduces the lattice correlator well on 
the largest lattice over rather wide ranges of bare quark 
masses and momenta. We use the new large volume re- 
sults to improve the extrapolation of the thermal quark 
mass values to the infinite volume limit. We obtain a 
value which is about 10% smaller in magnitude than the 
one obtained in the earlier analysis. The large spatial 
volume also allows one to directly analyze the momen- 
tum dependence of excitation spectra in more detail. We 
show that the dispersion relation for the plasmino mode 
obtained with the two-pole ansatz has a minimum at non- 
zero momentum. 

To clarify physical consequences for the quark spectral 
function that can be extracted from the analysis with 
the two-pole ansatz, we also take a closer look at the 
Euclidean quark correlator on the lattice. We show that 
with the present statistics the two-pole ansatz can give 
small x 2 /d°f even if the quark spectral function does 
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TABLE I. Simulation parameters on lattices. See the text 
for details. 

not consist solely of two peaks. Nevertheless, we argue 
that the analysis performed by us for finite ranges of 
parameters indicates the existence of quasi-particle peaks 
corresponding to the normal and plasmino modes having 
thermal masses near but above T c . 

At T = 1.5T C , we find that on about half of the gauge 
configurations the quark correlator shows an enhance- 
ment deviating from the generic behavior near the chiral 
limit. We show that this behavior is sensitive to the 
chirality of the quark propagator and argue that it is 
expected to arise from gauge configurations with non- 
trivial topology which are still abundant at this value of 
the temperature. We do not find any of these effects at 
the larger temperature T = 3T C . 

This paper is organized as follows. In the next section 
we summarize the basic properties of the quark corre- 
lator and present our simulation setup. In Sec. IIII1 we 
discuss the enhancement of the quark correlator near the 
chiral limit observed on some gauge configurations for 
T = 1.5T C . In Sec. II VI we analyze the dependences of the 
quark spectral function on the bare quark mass and mo- 
mentum. The last section is devoted to a short summary. 
In appendix [X] we discuss the relation between spectral 
function and Euclidean correlator to clarify the sensitiv- 
ity of the lattice correlator on the quark spectrum. 



II. QUARK SPECTRAL FUNCTION 

Excitation properties of the quark field are encoded in 
the quark spectral function p M „(u;,p), with /i and v de- 
noting Dirac indices. In order to extract p^iuj^p) from 
lattice QCD simulations we have analyzed the quark cor- 
relator in Euclidean space 

S^(r,p) = ^ J d 3 xd 3 y e*-<*-*><^(T, *)£,((), y)>, 

(1) 

on the lattice in the quenched approximation. This cor- 
relator is related to the spectral function as 

roo e (rT-l/2) u /r 

S^(r,p) = J dw cuj/2T - c _^ /2T PiiV (u,p). (2) 

Here, r is the imaginary time restricted to the interval 
< r < 1/T, ifj^TjX) is the quark operator, and V de- 
notes the volume of the system. The quark correlator 
Eq. (jTJ) has been calculated after fixing each gauge field 
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configuration to Landau gauge, d^A^ = 0. In the numer- 
ical analysis, we used (l/3)tr|9 M A M | 2 < 1CT 12 as a stop- 
ping criterion of the gauge-fixing algorithm. The simu- 
lation parameters are summarized in Table |TJ We used 
lattices of size N 3 x N T = 128 3 x 16 for T = 1.5T C and 
3T C . For the lattice fermion, we use non-perturbatively 
improved clover Wilson fermions with the clover coef- 
ficient csw taken from Rcf. 18 1. The number of con- 
figurations, iVconf, analyzed in the present study for 
T = 3T C is smaller than those in the previous work with 
-^conf = 44 — 60. With the aid of the definition of quark 
correlator with wall source [l8|, however, the statistics of 
the correlator improves compared to our previous analy- 
sis, as will be discussed in Sec. Hvl 

For T = 1.5T C , we find that the quark correlators on 
about half of the gauge configurations show an enhance- 
ment near the chiral limit deviating from the generic be- 
havior [2l[. As described in detail in the next section, 
this enhancement is obviously related to the non-trivial 
topology of gauge field configurations. The number of 
configurations which are identified to be topologically 
non-trivial is given in Table U as AT top . We found that the 
influence of topology on the quark propagator manifest 
itself in the quark correlator in scalar, pseudo-scalar, and 
tensor channels, whereas vector and axial- vector channels 
are not affected. Since all analyses in this study depend 
only on the quark correlator in the vector channel except 
for the analysis of the quark mass dependence given in 
Fig. [SJ which shows the result only for T = 3T C , we have 
performed these analyses using all gauge configurations 
without distinction of topologically trivial and non-trivial 
ones. For T = 3T C all configurations were identified to 
be topologically trivial. 

In the Wilson fermion formulation the quark mass, mo, 
is controlled by the hopping parameter, k. When we 
discuss the mo dependence of the quark spectrum in what 
follows, we use the relation 



mo = — log 
a 



1 If-- — 

2 \ K K r 



which defines the pole mass of the Wilson fermion propa- 
gator. In the free case, corresponding to ft — > oo, one has 
k c = 1/8. At finite values of the gauge coupling, (3 < oo, 
the quark mass receives an additive renormalization. In 
this work we determine k c from the behavior of the quark 
correlation function [l8[ , which is presented in Table fl] 



III. QUARK CORRELATOR ON 
TOPOLOGICALLY-NONTRIVIAL 
CONFIGURATIONS 

In the analysis of the quark correlator we found that 
for T = 1.5T C the correlator near the chiral limit shows 
a large enhancement deviating from the generic behavior 
fl8| on about half the gauge configurations analyzed in 
this study. At T = 3T C , on the other hand, no such effect 
has been observed. 
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FIG. 1. Real and imaginary parts of quark correlator 
in the pseudo-scalar channel, S-y 5 (n, T2), on a topologically- 
nontrivial gauge configuration for T = 1.5T C . 
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FIG. 2. Quark correlators in the temporal component of 
(3) vector channel, 5 7D (r), averaged over topologically trivial and 



non-trivial gauge configurations for T = 1.5T C 



In this section, we summarize the behavior of the quark 
correlator on these configurations. We show that char- 
acteristic features of this enhancement of the quark cor- 
relator are consistent with expectations for the behavior 
of a quark propagator on gauge field configurations with 
non-trivial topology which give rise to zero modes of the 
Dirac operator with definite chirality (22j . 

The Dirac structure of the quark correlator for p = is 
decomposed into its different quantum number channels 



Sr(ri,7s) = — 



d 3 xd 3 yTr D [r(^(n,x)V(7a,i/))] 



(4) 



where Ttq denotes the trace over Dirac indices with 
T = 1, 75, 7 M , 7^75, and i^, j„]/2 for scalar, pseudo- 
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total No N L Nr Nlr 
101 42 14 40 5 



TABLE II. Number of configuration for T = 1.5T C having no 
anomalous behavior (No), having in left- and right-handed 
channels (Nl and Nr, respectively), and in both channels 
(N LR ). 



scalar, vector, axial-vector, and tensor channels, respec- 
tively. In Fig. Q] we show the quark correlator in the 
pseudo-scalar channel, S* 75 (ti,T2), in the vicinity of the 
chiral limit on a configuration which is identified to be 
topologically non-trivial. The quark correlator in the 
pseudo-scalar channel should vanish when the system is 
invariant under parity transformation. In fact, on all con- 
figurations for T = 3T C , we have checked that S 75 (ti, t 2 ) 
is negligibly small compared to scalar and vector chan- 
nels. For T = 1.5T C , however, on about half of the con- 
figurations S* 75 (n, T2) shows an enhancement as in Fig.[T] 
near the chiral limit. The amplitude of this behavior 
rises with increasing k, i.e., decreasing quark mass, and 
takes the largest value for k near the critical value k c , 
at which the sign is flipped. A similar behavior is ob- 
served in scalar and tensor channels on these configura- 
tions. On the other hand, such an enhancement does 
not appear in vector and axial-vector channels. As an 
example, in Fig. [5] we show the temporal component of 
the quark correlators in the vector channel with p = 0, 
S~/ (t, 0), in the vicinity of the chiral limit averaged sep- 
arately over topologically trivial and non-trivial configu- 
rations for T = 1.5T C . The figure shows that the vector- 
channel correlator on each set of configurations does not 
give a statistically-significant difference. 

Our numerical result also shows that on many configu- 
rations which are identified to be topologically non-trivial 
the enhancement manifests itself only on either the left- 
or right-handed correlator, which arc defined as 



Sl,r(t,p) 



iTr D [A ± S(r,p)], 



(5) 



with projection operators onto left- and right-handed 
quarks, A± = (1 ± 7s)/2. In Table HU we show the num- 
ber of configurations having the enhancement in the left- 
and right-handed channels (Nl and Nr), respectively, 
and without the enhancement (No). On some configura- 
tions an enhancement of the correlator is observed both 
in left- and right-handed channels. Their number is de- 
noted by Nlr and is also given in Table |TTJ 

These results on the enhancement of the quark corre- 
lator near the chiral limit strongly indicate that they are 
related to a zero mode of the Dirac operator which, for 
instance, may be associated with the presence of instan- 
tons on these configurations. Using the eigenvalues and 
eigenmodes of the Dirac operator D = 7^(9^ + igA^), 



Dtpx(x) = Xip\(x), 



(6) 



the quark propagator, which is the inverse of iD + to, is 
written as 



S(x,y) = 



1 



iX + m 



(7) 



where x and y represent positions in Euclidean space. 
On a gauge field configuration having an instanton, the 
Dirac operator has a zero mode ipo [221 ] . and the quark 
propagator is decomposed as 



S(x,y) = — rlio(x)ip f (y) + Y] ^ — M^xiv) 

TO * — ' iX + TO 

A^O 

= Si(x,y) + S\(x,y), 



(8) 



with 



Si(x,y) = —ipo(x)ipl(y). 



(9) 



The zero mode associated with an instanton belongs to 
the left-handed spinor and satisfies A+Vo = t/j . Accord- 
ingly, Si satisfies 



Si(x,y) = A + Si(x,y)A A 



(10) 



On a configuration having an anti-instanton, the zero 
mode is right-handed, A-ipo = "00, and A + in Eq. (fT0|) 
are replaced with A_. Eq. (fTU)l shows that the zero mode 
cannot affect vector and axial-vector channels, since in 
these channels T in Eq. Q contains one gamma matrix, 
and hence Tr D [rA±] = 0. 

All these properties of Si agree with the behavior of the 
quark correlator observed on the lattice for T = 1.5T C 
discussed above. Gauge configurations in Nl and 
would have one or more instantons and anti-instantons, 
respectively, and those in Nlr have both an instanton 
and an anti-instanton. We have also checked that the 
quark correlator near the chiral limit on topologically 
non-trivial configurations approximately satisfies 



Sr(n,T 2 ) = [Sr(T2,n)\* 



(11) 



for pseudo-scalar and tensor channels (an example is 
given in Fig. [T]), which is consistent with Eq. 

Assuming that the enhancement of the quark correla- 
tor indeed arises from configurations with non-zero topo- 
logical charge one can estimate the topological suscepti- 
bility at that value of the temperature 



XQ 



T(SQ 2 
V 



(12) 



where Q = n+ — n- is the topological charge with n± 
being the numbers of instantons and anti-instantons, re- 
spectively, and VT -1 is the four- volume in Euclidean 
space. Provided that each configuration in A^ and Nr 
contains only one instanton and anti-instanton, respec- 
tively, and Q = on configurations in Nlr-, one obtains 
(SQ 2 ) = (N L + N R )/N conf ~ 0.53 for T = 1.5T C . Using 
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this result in Eq. (fT2"|) with VT 3 = 8 3 we obtain an es- 
timate for the topological susceptibility of Xq 4 — 0.1 IT 
for this temperature. We emphasize that this is a rough 
estimate since the value of (8Q 2 ) alters if some topo- 
logically non-trivial configurations have more than one 
(anti-)instanton. Moreover, while the separation between 
gauge configurations is a few times larger than the auto- 
correlation length measured in terms of the plaquette and 
Polyakov loop correlation functions [2l[ , this may not be 
the case for topology and the gauge configurations may 
not be separated well enough in Monte Carlo time to sup- 
press a long auto-correlation of Q. It is, however, notable 
that the value of \Q estimated here seems to agree with 
the previous study of \Q on quenched lattice for non-zero 
T [2|. 

The influence of non-trivial topological structures on 
the quark correlator is expected to be non- vanishing also 
in the thermodynamic limit on the quenched lattice. 
First, we have numerically checked that the effect does 
not cancel out on a single configuration when averaging 
correlation functions over the temporal direction, 



S(r) 



(13) 



Second, the numbers of (anti-)instantons, n±, increase 
proportional to V, while (|<5Q|) is proportional to I/ 1 / 2 
[24| . This insures that the topological susceptibility is 
well defined in the thermodynamic limit. Third, instan- 
tons and anti-instantons affect different channels in the 
quark correlator, respectively, and the numerical result 
indicates that their effect has a definite sign in each chan- 
nel. Therefore, effects of topology do not cancel out 
by their average in the thermodynamic limit. The non- 
trivial topological properties of the quark propagator ob- 
served here thus will give rise to physical contributions 
in various observables, such as those observed as U(1)a 
anomaly. 

Although the non-trivial topology of gauge configu- 
rations gives rise to non-vanishing contributions to the 
quark correlator, our numerical result for T = 1.5T C 
shows that the quark propagators in the vector and axial- 
vector channels are not affected by this effect as shown in 
Fig. O Since all analyses in Sec. II VI and appendix IA1 rely 
only on the vector channel, except for the one for Fig. [SJ 
we use all gauge configurations including topologically 
trivial and non-trivial ones in these analyses. ^From the 
insensitivity of the vector and axial-vector channel cor- 
relators on topology, it is anticipated that the instan- 
ton distribution for T = 1.5T C is sufficiently dilute so 
that each instanton affects the quark propagator as Si 
in Eq. ([5]) almost independently |24[. The analysis on 
the thermal mass and the dispersion relation performed 
in Sec. IIVI would then be interpreted as the one for the 
remaining part S\ . 

The discussion on the relation between the topology of 
gauge configurations and the enhancement of the quark 
correlator given in this section clearly addresses an im- 



portant aspect of QCD thermodynamics and requires fur- 
ther analysis. We intend to explore the interplay of topol- 
ogy, zero modes, eigenfunctions and the structure of the 
quark propagator in a forthcoming publication (25j . 



IV. QUARK EXCITATION SPECTRA 
A. Quark spectral function and two-pole ansatz 

The Dirac structure of the quark spectral function 
p^ v {ui,p) at finite temperature is decomposed as 

= p (w,_p)(7°) AI „ - p v (uj,p)(p ■ j) MV + p a (w,p)ln V , 

(14) 

where p = \p\ and p = p/p [2|. In this section we consider 
the spectral function above T c for two cases; (1) for zero 
momentum, and (2) in the chiral limit fl8j . With p = 
0, p v (uj,p) vanishes in Eq. (fT4|) and p^(uj,p = 0) is 
decomposed with the projection operators L± = (1 ± 
7°)/2 as 



p(w, 0) = pf ( W )L +7 U + ^M£-7 U , (15) 



with 



p^H = -Tr D [p( W> 0)7°I±] 



(16) 



In the chiral limit and for T > T c , the system possesses 
the chiral symmetry and p s (uj,p) vanishes. p^ v {ui,p) is 
then decomposed with the projection operators P± (p) = 
(l±7°p-7)/2 as 

p(u,p) =pP( W ,p)P + (p) 7 ° + pP( W ,p)P_(p)7 , (17) 

with 

pi (oj, P ) = iTr D [ P (cj, p)j°P ± (p)] . (18) 

Using the charge conjugation symmetry one can show 
that 

p¥(u) = p™(-u), pl^p) = pP(- W ,p). (19) 

In the chiral limit and for p = 0, both p s (uj,p) and 
p v (uj,p) vanish and one obtains 



p£M=/4(w,o) = po(«,o). 



(20) 



In order to extract the quark spectral function from 
lattice correlator, we follow the approach taken in 
Refs. [13, [3: which makes use of a two-pole ansatz for 
the spectrum, 



pf P (lu) = ZiSiu - E x ) + Z 2 5{u + E 2 ). (21) 
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FIG. 3. Lattice correlator So(T)iat in the chiral limit nor- 
malized by the fitting function Eq. (gU for T/T c = 3. Solid 
and dashed lines represent the results for N a /N T — 8 and 4, 
respectively. 
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FIG. 4. T m i n dependence of fitting parameter E in the fit 
with Eq. ((22)1 to Sq(t, 0)i at in the chiral limit. 



Here, ^1,2, and Ex^ > 0, are fitting parameters that 
will be determined from correlated fits to the lattice cor- 
relator: Zi t 2 and E\ y 2 represent the residues and posi- 
tions of poles, respectively. By comparing the fit results 
with spectral functions obtained in perturbative calcu- 
lations one can identify the pole at uj = E\ to be the 
normal mode, while the one at w = — E2 corresponds to 
the plasmino mode @, H, [H| ■ To determine the fit pa- 
rameters with correlated fits, we use lattice data points 
at r m i n < t /a < N T — T m in with the number of points 
iVdata = N T - 2r min + 1. 

Before discussing the analysis of the quark spectrum 
with this fitting function, let us examine the behavior 
of the quark correlator, Eq. JT]), on the lattice and dis- 
cuss how the fitting function Eq. (j2~Tj) reproduces it. To 
simplify the argument we limit our attention for the mo- 
ment to the case in the chiral limit with p = 0. The Dirac 
structure of the spectral function is then proportional to 



7 , p(uj,0) = ^0(^,0)7°, and the Euclidean correlator 
is solely given by the one corresponding to this channel, 
So(t, 0) = Ttd[S(t, 0)70]. Since Po(uj, 0) is even in lu, the 
ansatz Eq. (|^T|) reduces to 



Po(^,0) 



Z(6(u 



~ P >,0) 



E) + S(lo + E)) , (22) 



which involves only two fit parameters, Z and E. In 
Fig.[3jwe show the lattice correlator Sq(t, 0)i a t for T/T c = 
3 normalized by the one constructed from the fitting 
function Eq. (|2"2"]) . Sq(t, 0)at, with fitting parameters de- 
termined from a correlated fit with r m - m = 4. The fig- 
ure shows that the fit result, Sq(t, 0)fit, well reproduces 
the lattice results, Sq(t, 0)i at , for tT ~ 0.5. We have 
checked that the fitting parameters do not have a sta- 
tistically significant dependence on the choice of r m i n for 
3 < Tmin < 6. This is evident from the dependence of E 
on T m ; n . Results obtained in the chiral limit are shown 
in Fig. SI In Fig. [31 clear deviations between Sq(t, 0)at 
and Sq(t, 0)i a t are observed for small and large r, which 
can be attributed to distortion effects arising from the 
source [H, Hfl- For r min < 2, x 2 /dof of the fit becomes 
unacceptably large due to the distortion effect. Similar 
results are obtained for all parameters analyzed in this 
study with Eq. (|2f p . In the following, we use T m - m = 4, 
which is chosen in order to reduce the distortion effect 
and at the same time leave us with a sufficient number 
of data points for our four parameter fits. We found that 
with r m j n = 4 fits with the pole ansatz give reasonable 
values of chi-square, 0.5 < x 2 /dof < 3, for m /T < 0.5 
and p/T < 3. 

We remark that, although the two-pole ansatz Eq. (|2Tj) 
reproduces the lattice correlator quite well with a small 
X 2 /dof, this result does not necessarily mean that the 
quark spectral function is solely composed of two sharp 
peaks at energies u = E\ and — E%. In fact, as discussed 
in appendix \K\ the Euclidean correlator is insensitive to 
the structure of the spectral function at \u>\ < T. Al- 
though with our current statistics the quark propagator 
at all distances is determined to better than 0.1%, the 
spectrum in this low energy range generally still suffers 
from large uncertainties. In Sec. IIV CI and appendix \^ 
we will discuss how our analysis of the quark spectrum is 
influenced by this problem. Nevertheless, we will assert 
in these sections that, from the success of the two-pole 
ansatz over wide ranges of mo and p, it is reasonable to 
expect that the quark spectrum contains peak structures 
corresponding to the normal and plasmino modes. Qual- 
itative dependences of the positions of these peaks on Too 
and p are then analyzed with the fitting parameters E\ 
and E 2 - 

Finally, we briefly comment on the statistical errors 
of the quark correlator in this analysis. In Fig. [3J the 
quark correlator obtained in Ref. [18| with N a /N T = 4 
for T = 3T C is shown by the dotted lines together with 
the present result with N a /N T = 8. The figure shows 
that the error-bars for N a /N T = 8 is about 40% smaller 
in magnitude than that for N a /N T = 4, although the 
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FIG. 5. Bare quark mass dependence of parameters Ei, Ei 
at T = 3T C on lattices with N a /N T = 8, 4 and 3. 



N a /N T 


3 4 8 oo 


T/T c = 3 
T/T c = 1.5 


1.041(16) 0.875(8) 0.743(6) 0.725(14) 
1.075(14) 0.906(8) 0.784(3) 0.768(11) 



TABLE III. The values of ttit/T obtained on lattices with 
N a /N T = 3, 4 and 8. The far right column shows the values 
of vtlt estimated by the extrapolation to the infinite volume 
limit. 



former result is obtained with only about half as many 
configurations than the latter (iV con f = 51). As men- 
tioned previously, this improvement in statistics is a con- 
sequence of the use of the wall source in the numerical 
analysis of quark correlators |l8[ . We have checked that 
the distribution of correlators obtained on each gauge 
configuration around the average is approximately pro- 
portional to V^ 1 ! 2 . Such a scaling with the wall source 
is to be expected when the spatial extent of the lattice, 
L, is larger than the coherence length of the quark field, 
£. Estimating £T ~ T/m T ~ 1.25 and LT = N a /N T = 8 
this, indeed, seems to hold. 



B. Excitation spectra for p = 

In this subsection, we focus on the quark spectrum 
with p = 0, which is decomposed into p±{u>) as in 
Eq. I|15p. and analyze the mo dependence of the spec- 



trum. In Fig. [5l we show results for the too dependence of 
fitting parameters, E\, E2, and Z2I (Zx + Z^) for T = 3T C . 
We compare the results on lattices with smaller aspect 
ratios N a /N T = 3 (triangles) and 4 (squares) [UEII with 
our new results on the 128 3 x 16 lattice (circles). One sees 
that the qualitative behavior of fitting parameters on the 
largest lattice agrees with previous results obtained on 
lattices with smaller N a /N T . The ratio Z^jiZx + Z2) be- 
comes larger as m is decreased, and eventually reaches 
Z2KZ1 + Z%) = 0.5 in the chiral limit. At this point, 
which defines k c in Eq. ([3]), E\ = E<i is satisfied within 
the statistical error. The thermal mass is defined as 
rriT = (Ei +Ez)/2 on this point, or equivalently ttit = E 
with the fit with Eq. (|25J) to 5 (r, 0). For m > 0, E x has 
a minimum, while E2 is a monotonically increasing func- 
tion. For T = 1.5T C , the same analysis can be performed 
on topologically trivial gauge configurations, which gives 
qualitatively the same result as in Fig. [5j 

Figure [S] also shows that the values of E\ and E2 drop 
significantly with increasing aspect ratio N a /N T . Ac- 
cordingly, the value of mj- also has a strong N a /N T de- 
pendence. The value of ttit/T for N a /N r = 3, 4 and 8 
for T/T c = 1.5 and 3 is given in Table [TTH To infer the 
thermal mass in the infinite volume limit, wc performed 
an extrapolation of mj to infinite volume with an ansatz 



mr(l/V) - m T (0)exp(c/V) 



(23) 



using the results with N a /N T = 4 and 8. The result of 
this extrapolation is shown in the far right column of Ta- 
bic IIIII The infinite volume extrapolated values of tut 
for each T obtained in this study is about 10% smaller 
in magnitude than those estimated in Rcf. [18]. It is 
also notable that the ratio ttit/T increases as T is low- 
ered, which was not observed in the previous study [l8|. 
With the present extrapolation, the value of tot in the 
infinite volume limit coincides with the one obtained for 
N a /N T = 8 within the statistical error. This result sug- 
gests that the influence of finite volume effects on itlt is 
well suppressed for N a /N T = 8. 

In Fig. one also finds that the value of k c , i.e. Too 
at which Z2/(Z\ + Z2) = 0.5, slightly shifts with the 
variation of N a /N T . The difference of k c between the 
two largest lattices N a /N T = 4 and 8, however, is within 
the statistical error. 



C. Dispersion relations at non-zero p 

Next, we set k = k c and analyze the momentum depen- 
dence of the excitation spectra at non-zero momentum 
in the chiral limit using the decomposition in Eq. (|17p . 
For the analysis of the quark correlator at non-zero p for 
T = 1.5T C , we have performed the analysis on 54 con- 
figurations including topologically trivial and non-trivial 
ones. In Fig.[6]we show the momentum dependence of E\ 
and E2 normalized by tot, as well as ^2/(^1 + ^2), for 
T/T c = 1.5 and 3. The horizontal axis represents the mo- 
mentum of a free fermion on the lattice, p = (1/a) sin pa, 
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FIG. 6. Dependences of the fitting parameters E\ and E2 
and the ratio Z2KZ-1 + Z2) on the lattice momentum p — 
(l/a)sin(pa) for T/T c — 1.5 and 3. See, the text for the 
explanation of other lines. 



normalized by mx- The figure shows that the qualita- 
tive features of the dispersion relations as functions of p 
do not change compared to our earlier results on lattices 
with smaller aspect ratio [l8| : Ei > E 2 is always satisfied 
in contrast with the results in the previous subsection, 
and E2 enters the space-like region at high momentum. 
We note that a similar result is obtained from the analysis 
of the quark spectrum in the Schwinger-Dyson approach 
[IH, [l3|. The decrease of Z 2 j(Z\ + Z 2 ) for large momen- 
tum shows that the plasmino mode ceases to exist as p 
becomes larger. 

At asymptotically high temperatures the quark spec- 
trum is calculated perturbatively @, [2(| ■ The dispersion 
relations of the normal and plasmino modes obtained at 
the one- loop order are shown in Fig. [6] by the dotted lines 
with rriT = gT/^/6. The residue of the plasmino mode in 
the normalization of the quark field /_ dujp±(u),p) = 1 
is also shown in the lower panel by the dotted line. Be- 
sides these poles, p+(u;,p) at the one-loop order has a 
continuum in the space-like region. The spectral weight 
of the continuum, Z cont . = duj p^_(uj , p) , is shown in 
the lower panel by the dashed line, which in the pertur- 
bative analysis reaches at most exceeds Z COTlt ~ 0.18. 

The numerical result for E2 presented in Fig. [5] shows 
that the result for the lowest non-zero momentum, 
Pmin = 27r(iV T /iV (T )T, is significantly lower than m^; 
i. e. m T /T ~ E 2 /T\ p=0 = 0.742(12) and E 2 /T\ p=Pmixi = 
0.694(21) for T/T c = 3 and E 2 /T\ p=Q = 0.782(9) and 



E 2 /T\ p = Pmia = 0.758(18) for T/T c = 1.5. Provided that 
the value of E 2 in our two-pole ansatz represents the dis- 
persion relation of the plasmino mode, this result serves 
as direct evidence for the existence of the plasmino min- 
imum in the non-perturbative analysis. One, however, 
has to be careful with this interpretation. As we have 
argued previously and in appendix \^ the Euclidean cor- 
relator is insensitive to the spectral function at low en- 
ergy, \uj\ < T, and analysis of the spectrum in this energy 
range has large uncertainty. In fact, in terms of the com- 
parison with the perturbative result our fitting function 
does not contain the continuum in the space-like region. 

In order to see the effect of the continuum part in 
the spectral function at the one-loop order on the Eu- 
clidean correlator, we have performed the following anal- 
ysis. First, we calculated the Euclidean correlator at 
one- loop order in perturbation theory, S*htl (t, p) . We 
then performed the two-pole fit with Eq. (|2"Tj) to discrete 
points of ShtiX 7 ", p) with the same discretization as on 
the lattice. To estimate the chi-square of the fit, we put 
error-bars and Gaussian noise on each value of 5*htl (t, p) 
following the covariance matrix obtained on the lattice. 
With this analysis, we found that the two-pole fit to 
Srtl(t,p) gives a reasonable chi-square, x 2 /dof ~ 1, al- 
though the spectral function corresponding to S*htl (t, p) 
contains the continuum. The values of E\ and E 2 ob- 
tained by the fit to 5 , htl(t, P) are shown by the dashed 
lines in the upper panel of Fig. [6] One sees that these 
results considerably deviate from the positions of normal 
and plasmino modes. In particular, we note that E 2 also 
enters the space-like region for p > T. This result obvi- 
ously shows that the success of the pole ansatz does not 
mean the absence of the continuum. The values of E\ 
and E 2 thus should not be regarded as the positions of 
peaks, but a center of spectral weights in some sense in 
the positive and negative energy regions. It also suggests 
that in our current analysis the normal mode is quan- 
titatively much better under control than the plasmino 
branch. 

Nonetheless, some qualitative features on the quark 
spectrum can be extracted from our results. First, as 
discussed in Appendix [2] the consistent success of the 
two-pole ansatz over wide parameter ranges suggests that 
the quark spectrum has peak structures at positive and 
negative energies in the range of mo and p analyzed in 
this study. Next, the momentum dependences of the 
peak positions in fPAui^p) can be constrained by gen- 
eral properties of the spectral function. Using the rela- 
tion Eq. (fT9l and the analyticity of p(ur,p), one can show 
that (1) if p^(u!,p) at p = has a peak at positive energy 
uj = E+(0), there also exists a peak at negative energy 
with the same absolute value, oj = — E*_ (0) = — £+(0), 
and (2) momentum dependences of the peak positions 
uj = EUp) and -E*_{p) satisfy dE* + /dp = -dE*_/dp at 
p = 00. The result in Fig. [6] indicates not only a nega- 



A discussion of a similar relation for position of poles in the 
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tive dEz /dp at p = but also a positive dE\ /dp. It thus 
seems reasonable to conclude that the dispersion relation 
of the plasmino peak, E*L(p), has a negative derivative 
near p = 0. This result then leads to the existence of a 
minimum of E*_ (p), provided that the plasmino peak sur- 
vives at sufficiently high momentum where the dE*_/dp 
becomes positive. 



V. BRIEF SUMMARY 

In this study, we extended the analysis of the quark 
spectral function on quenched lattices, performed in pre- 
vious work [It], [H!]j to lattices with larger spatial volume. 
With this analysis we achieved considerable progress in 
the understanding of finite spatial volume effects on the 
quark correlator and confirmed the importance of non- 
trivial topological configurations for the structure of the 
quark propagator at temperatures above T c . The strong 
V dependence observed in a previous study with aspect 
ratios N a /N T = 3 and 4 seems to converge on the largest 
lattice with N a /N T = 8. An extrapolation of the quark 
thermal mass, mr, to the infinite volume is about 10% 
smaller than the earlier estimate [l8[ . 

Dependences of the fitting parameters in the two-pole 
ansatz on mo and p do not change qualitatively on the 
largest lattice. We, however, found that the position of 
the pole corresponding to the plasmino, E2, has a mini- 
mum at non-zero momentum within the ansatz. As dis- 
cussed in Sec. IIV C[ this result with the pole ansatz does 
not directly mean the existence of the minimum of the 
plasmino dispersion (See also, Ref. [Hf-)- The quark cor- 
relator obtained on the lattice as functions of continuous 
parameters mo and p, however, strongly indicates the 
existence of the plasmino peaks and that the dispersion 
relation of the plasmino has a negative slope at low mo- 
mentum near T c . 

We also presented evidence for the influence of non- 
trivial topology of gauge field configurations on the struc- 
ture of the quark propagator. This clearly is visible in our 
data at T = 1.5T C , while at T = 3T C we do not observe 
topologically non-trivial configurations. We showed that, 
nonetheless, the vector and axial vector channels are not 
influenced by this effect at T = 1.5T C . We expect that the 
effect of non-trivial topology will modify the structure of 
the quark propagator in the quenched QCD significantly 
as T approaches T c from above. The analysis of correla- 
tors near T c therefore requires a careful treatment of this 
effect. 



quark spectral function is presented in [27f The extension of this 
argument to the position of peaks is straightforward. This feature 
is graphically seen, for example, in the contour maps of spectral 
functions in Ref. 



Appendix A: Difficulty in analyzing the low energy 
structure of spectral functions 

In this study, we have used the two-pole ansatz Eq. (|2Tj) 
to examine the quark spectrum from Euclidean lattice 
correlator. This ansatz is found to reproduce the lattice 
correlator well over wide parameter ranges. In this ap- 
pendix, we address the significance of the success of the 
two-pole ansatz. We show that this success to some ex- 
tent is due to a poor resolution of the lattice correlator on 
the spectral function at low energy, |w| < T. The insensi- 
tivity of the Euclidean correlator to low-energy spectrum 
has been recognized through the studies of spectral func- 
tions on the lattice, especially efforts to extract transport 
properties psT - lsoj . Since the spectral weight of quarks 
with a small mass and a momentum concentrates in the 
energy range \u\ < T, our analysis suffers from the same 
difficulty. 

In the next subsection we first summarize the problem 
of the insensitivity in terms of the power series expansion 
of correlators. We then take a closer look at the quark 
correlator in Sec. IA 21 

1. Moment expansion 

Let us consider a Euclidean correlator S (r) and a spec- 
tral function p{ui) which are related to each other as 

roo e (-rT- 1/2 V/T 

S{T) = J-00 ^ e^ + C^ 2T P ^' (M) 

where £ = ±1 for fermions and bosons, respectively. By 
Taylor expanding exp[(rT— l/2)w/T], Eq. ([AT) is written 
as 

S{r)=Y,C n {rT- l X , (A2) 

where the coefficients C n are given by the moments 
of thermal spectral function, p'(w) = p(uj)/(e ul ^ 2T + 

Cn = - J <*"(y) P'(",P)- (A3) 

One can show that the power series Eq. (|A2|) converges 
for < t < 1/T for spectral functions having an asymp- 
totic form p(uj) ~ (uj rn for large |w|. 

In lattice simulations one obtains the values of S(r) 
for N T discrete r values with statistical errors. Through 
Eq. (|A1[) . these values of S(t) provide N T 'pieces of dif- 
ferent' information on the structure of p(uS). The coeffi- 
cients of the power series, C„, can be interpreted as the 
different representation of this information. 

To consider how the information of low-energy struc- 
ture of p(ui) is encoded in S(t), we decompose p(ui) into 
low and high energy parts as p(oj) = pa(lo) + Phigh(w) 
with pa(w) = p(u))6{h — \u)\) with some energy A. Since 
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S(t) and p(to) arc related with each other linearly, one 
can then also decompose S(t) into contributions of low 
and high energy spectrum. The contribution of Pa(w) is 
written as 



s a{t) = / dw 



e (rT-l/2) W /T 
j/2T _|_ ^ e -w/2T 



with 



u.' 



A 

-A VT 



//(a;). 



(A4) 

(A5) 
(A6) 



To avoid unnecessary complexity in the following we 
limit our attention to the spectral function satisfying 
p(u>) = Cp(-u>) > for cj > 0. Odd terms in Eq. (|A6]) 
then vanish while C4 > for even n. With the aid of in- 
equalities (w/T)" < (A/T)" in the integrand of Eq. (TA~6|) 
and |tT — 1/2| < 1/2, one then obtains an inequality 



sL A \r) Sr(r) 1 

" ra! V 2T 



S(r) 



A 



(A7) 



This inequality provides an upper limit of the contribu- 
tion of Sn (t) to S'(t), ie. the relative strength of sig- 
nals for information on each moment of p'^{oj) encoded in 
S(t). Eq. (|A7j) tells us that if we choose A < T the con- 
tribution of higher order moments is strongly suppressed 
due to the factorial and exponential factors. The decrease 
of the signals makes evaluations of their values difficult. 
In practical analyses with statistical errors, therefore, the 
number of moments, i.e. the number of independent in- 
formation on pa(w), which can be estimated from the 
correlator is severely limited with Eq. (|A7[) . For exam- 
ple, if one wants to examine /Oa(w) with A = T from the 
lattice correlator having statistical errors of order 0.01%, 
S { 6 A \t)/S(t) < (l/6!)2- 6 ~ 2 x 10~ 5 is already signifi- 
cantly smaller than the statistical error. Moments of or- 
der higher than the sixth are inaccessible unless the cor- 
relators with different r values have strong correlations 
corresponding to the moments. One can obtain at most 
three statistically meaningful information on p\(uj), i.e. 

Cq A \ C*2 A ' and C± , from the correlator. We empha- 
size that this limitation on the number of independent 
information on pa(w) is not resolved with the increase of 
N T . 

When phigh(w) gives rise to a large contribution to 
S(t), the analysis of pa(w) should be much more diffi- 
cult. This occurs particularly when p{uj) has a positive 
mass dimension, m, since p(w) behaves as ~ w m for large 
ui for this case. Since the quark spectral function has a 
negative mass dimension, effects of /9high(w) is expected 
to be small in our analysis. 

While in the above argument we have limited our at- 
tention to the spectrum satisfying p{uj) = C,p{—io) 1 one 
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FIG. 7. Mid-point subtracted correlator So(t, 0) in the chi- 
ral limit (upper panel), and So(t,p)/(tT — 1/2) 2 for p = 
(middle) and 2(2nN T /N a )T (lower). 



can also derive a similar inequality for odd terms of 
S^\t) for general cases. 



2. Quark correlator 

Next, we inspect the structure of the quark correlator 
based on the discussion in the previous subsection. Let 
us first consider the correlator in the chiral limit with 
p = 0. In this case, the quark correlator takes non- 
zero values only in the temporal component of the vec- 
tor channel, Sq(t, 0)i at . The correlator in this channel 
is symmetric, Sq(t, 0) = Sq(1/T — r, 0), and odd terms 
vanish in Eq. (|A2I) . The values of C n for even n are eval- 
uated from the correlator as follows. First, the lowest 
moment, Co, is directly read off from the correlator at 
t = 1/(2T); Co = So(l/(2T),0)i at . To examine higher 
order moments, we show in Fig. [7] the r dependence of 
the mid-pofjpfcGsubtracted correlator 



S (t,0) = S (r,0) - S (1/2T,0), 



(A8) 



in the upper panel and that divided by (tT — 1/2) 2 
in the middle panel. Using the latter plot, C2 is esti- 
mated as C 2 = lim rT _,. 1/2 5o(T,0)/(rr - 1/2) 2 , while 
the deviation from the constant in this plot carries in- 
formation of C„ for n > 4. The figure, however, shows 
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that Sq(t,0)/(tT — 1/2) 2 is constant within the statis- 
tical error except for the range of r near the source, 
\tT — 1/2 1 > 0.3, which receives the distortion effect 
and should be excluded from the analysis. This result 
shows that information on moments higher than four are 
buried in the statistical errors and inaccessible from our 
numerical results for this parameter. 

This result is consistent with the argument in the pre- 
vious subsection. As the pole ansatz predicts the poles 
of po(w,0) at |w|/T ~ 0.74, spectral weight of po(uj,0) 
is expected to concentrate more or less in the energy 
range \uj\ < T. Substituting A = T to Eq. (|AT|) . one 

obtains that S^\t, 0)\&t/ S(t, 0)i a t arc suppressed faster 
than (2"n!) _1 . This upper limit is already the same order 
as the typical magnitude of the statistical error of the cor- 
relator, 10~ 3 , for n = 4. One thus can obtain only two 
statistically meaningful information on the quark spec- 
trum, Co and C2, from the correlator for this parameter. 

Since the quark correlator for mo = p = contains 
only two information on the quark spectrum, any ansatz 
having more than two fitting parameters can well repro- 
duce the correlator as long as it gives sufficiently small 
higher order moments. For example, we found that a fit- 
ting function with a single Gaussian peak for pq(lo,0) in 
the chiral limit, 

po(w,0) = #M = -|-cxp f-^j , (A9) 

including two fitting parameters, Z and T, gives 
X 2 /dof ~ 2.5 for T/T c = 3, which is less than twice 
larger than the one with the two-pole ansatz. 

Evaluation of the quark spectrum for too = p = 
with the lattice correlator, therefore, is quite difficult. 
It, however, is notable that one can obtain some addi- 
tional information on the spectrum for too = p = from 
the numerical analysis by using correlator for non-zero 
too and p, because the spectral function is continuous 
with respect to these parameters. When too (p) takes 
non-zero values, (p^_(uj,p)) is no longer an even 

function. Odd moments in Eq. (| A2|) then take non-zero 
values and these values provide additional constraints for 
the spectrum. At the same time the strength of /9+' P 
shifts toward higher energie s, as is expected from the 
free quark spectrum, E = \Jp 2 + tuq. The higher order 



moments can then become larger, which makes the esti- 
mates of their values easier. To see this explicitly, in the 
lower panel of Fig. Owe show So{r,p)\ a t / {tT —1/2) 2 with 
p = 2p min = 2(2nN T /N !y )T. The figure shows that this 
function clearly deviates from a constant, and hence con- 
tains statistically significant information on the fourth 
moment and higher. The number of information acces- 
sible with the quark correlator with this parameter thus 
is at least more than five. Since this number is larger 
than that of fitting parameters in the two-pole ansatz, the 
small x 2 /dof with the ansatz for this correlator is non- 
trivial. Since the simple two-poles ansatz Eq. (|2ip can 
reproduce the lattice correlator not only for too = p = 
but also over wide ranges of continuous parameters, it is 
expected that this ansatz reproduces at least a qualita- 
tive feature of the global structure of the quark spectrum 
in this parameter range, including mo = p = 0. 

Similarly, one can argue that the single Gaussian 
ansatz, Eq. (|A9[) . in the chiral limit and for p = is 
most probably excluded as follows. The continuity of 
the quark spectrum means that if the Gaussian ansatz 
appropriately reproduces the structure of the spectrum, 
the ansatz has to be continuously connected to non-zero 
mo. We, however, have checked that several ansatze for 
Pt:(co) which reduce to Eq. (|A9[) at mo = cannot give 
a reasonable chi-square with mo 7^ 0; as far as we have 
checked, x 2 /dof obtained on all the ansatze grows rapidly 
as mo is increased. Although this analysis cannot com- 
pletely exclude the form of Eq. (|A9[) in the chiral limit, it 
strongly suggests that the quark spectrum like Eq. (|A9[) 
is most probably ruled out even in the chiral limit. 
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